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1.0  INTRODUCTION 


The  ballistic  wind  error  analysis  models,  RAWIN,  RADAR,  and  NAVAID,  are  for¬ 
mulated  in  Volume  I  of  this  report.  The  reader  is  referred  to  Volume  I  for 
general  discussion  of  each  model  and  for  definitions  of  various  terms,  such 
ballistic  line,  component  velocity  variance  (CW) ,  etc. 

Each  error  analysis  model  has  been  programmed  in  the  ASCII  FORTRAN  language. 
Each  separate  program  is  named  for  its  corresponding  model.  The  implementa¬ 
tion  of  these  programs  on  the  UNIVAC  1108  computer  at  White  Sands  Missile 
Range  is  described  in  this  Users'  Manual. 

The  utilization  of  the  auxiliary  program  LRDC  is  described  in  Section  10. 


2.0  PRELIMINARY  CONSIDERATIONS 


Each  of  the  error  analysis  programs  RAWIN,  RADAR,  and  NAVAID,  consists  of  a 
single  main  program  unit.  No  subroutines  are  called.  Each  program  may  be 
executed  in  either  batch  or  demand  mode.  The  reader  is  referred  to  the  accom¬ 
panying  listings  for  actual  program  source  code. 

Execution  is  similar  for  all  three  programs.  One  complete  execution  of  any  of 
the  programs  is  referred  to  here  as  a  complete  run.  A  complete  run  consists 
of  one  or  more  separate  problems.  The  general  sequence  of  input  and  execution 
is  described  below. 


STEP  EXPLANATION 

1  Three  card  images  are  read  on  logical  unit  105.  These  are  used  to 
document  the  complete  run. 

2  Each  disk  file  described  in  Section  4  is  read  once  on  its  appro¬ 
priate  logical  unit.  The  program  uses  these  inputs  to  perform 
preliminary  computations  required  for  the  complete  run. 

3  Two  card  images  are  read  on  unit  105  in  order  to  initiate  the  first 
problem.  The  first  card  image  documents  the  problem,  while  the 
second  contains  data.  Execution  then  procedes.  On  .put  is  written 
to  logical  units  106  and  1020. 

4  Step  3  may  be  repeated  for  any  number  of  problems.  Program  execu¬ 
tion  is  terminated  only  by  substituting  an  end-of-file  image,  i.e., 
@E0F,  for  any  of  the  cards  in  Step  1  or  3. 

None  of  the  cards  noted  above  may  be  omitted.  In  demand  usage  each  of  the 
input  card  images  is  solicited  by  the  executing  program. 


None  of  the  programs  requires  more  than  6000  words  of  storage  (IBANK  plus 
DBANK)  on  the  UNIVAC  1108.  The  demand  time  required  for  the  execution  of  a 
complete  run  consisting  of  a  single  problem  is  on  the  order  of  ten  seconds. 


3.0  CARD  INPUT 


The  card  image  inputs  described  below  for  each  program  are  entered  on  logical 
unit  105.  Currently  105  is  taken  to  be  unit  5,  which  on  the  UN  I  VAC  1108 
corresponds  to  a  card  reader  in  batch  usage  or  to  terminal  keyins  in  demand 
usage.  Logical  unit  105  can  be  respecified  to  a  different  value  by  changing 
the  appropriate  data  statement  in  each  program. 


3.1  RAWIN  (Logical  Unit  105) 


CARD 

DATA 

FORMAT 

1 

COMM 

(20A4) 

COMM  contains  up  to  80  characters  which  are  used  to 
document  the  complete  run.  A  blank  card  may  be  used 
if  no  documentation  is  desired. 

2 

COMM 

(20A4) 

CO  MM  contains  up  to  80  characters.  It  may  contain  a 
user  supplied  list  of  input  files  used  for  the  com¬ 
plete  run  and/or  it  may  contain  further  comment. 
This  card  may  also  be  blank  if  desired. 

3 

INTR 

(A4) 

INTR  may  have  the  value  YES  or  NO  and  must  begin  in 
column  1.  If  INTR  is  YES,  results  of  preliminary  and 
intermediate  computations  are  outputted  to  unit  1020. 

If  INTR  is  NO,  only  final  results  are  outputted  to 
unit  1020. 

4 

COMM 

(20A4) 

COMM  contains  up  to  80  characters  which  are  used  to 
document  the  first  problem.  This  card  may  be  blank 
if  desired. 

5 

BEL,  RE,  BA,  RA,  RLED ,  RLEA,  FE 

Free  Field 

This  set  of  nonnegative  real  variables  is  input  for 
the  first  problem. 


BEL  is  the  bias  error  (degrees)  in  elevation 

tracking  of  the  apparent  target. 

RE  is  the  random  error  (degrees)  in  elevation 

tracking  of  the  apparent  target. 

BA  is  the  bias  error  (degrees)  in  azimuth 

tracking  of  the  apparent  target. 

RA  is  the  random  error  (degrees)  in  azimuth 

tracking  of  the  apparent  target. 

RLED  is  the  random  error  (meters)  associated 

with  the  measurement  of  the  displacement  of 
the  launch  site  of  the  balloon-radiosonde 
combination  from  the  receiving  set. 

RLEA  is  the  random  error  (degrees)  associated 

with  the  launch  azimuth. 

FE  is  the  foreground  elevation  (degrees). 

With  changed  information  or  values,  cards  4  and  5  may  be  repeated  in  sequence 
any  number  of  times  in  order  to  execute  further  problems.  Values  of  all 
variables  listed  for  card  5  must  be  entered  for  each  problem,  even  if,  for 
example,  only  one  value  is  changed. 


To  terminate  execution,  enter  a  final  input  card  containing  the  end-of-file 
image  @EOF.  If  this  is  not  done,  the  program  will  expect  comment  and  data 
pertaining  to  further  problems. 


3 . 2  RADAR  (Logical  Unit  105) 

CARD  DATA  FORMAT 

1,2,3  Same  description  and  formats  as  for  RAWIN 

4  Initiates  first  problem.  Same  description  and  for¬ 
mats  as  for  RAWIN. 

5  BEL,  RE,  BA,  RA,  BS,  RS,  RLED,  RLEA,  FE  Free  Field 


This  set  of  nonnegative  real  variables  is  input  for 
the  first  problem. 
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BEL,  R£ ,  BA,  RA,  RLED,  RLEA,  and  FE  are  the  same  as 
described  for  RAWIN. 


BS  is  the  bias  error  (meters)  in  slant  range. 

RS  is  the  random  error  (meters)  in  slant  range. 

Further  problems  are  initiated  in  the  same  manner  as  described  for  RAWIN. 

Input  of  the  card  image  @EOF  is  required  to  terminate  execution. 

3 . 3  NAVAID  (Logical  Unit  105) 

CARD  DATA  FORMAT 

1.2.3  Same  description  and  formats  as  for  RAWIN, 

A  Initiates  first  problem.  Same  description  and 

mats  as  for  RAWIN. 

5  REX,  REY,  RLE  Free  Field 

This  set  of  nonnegative  real  variables  is  input  for 

the  first  problem. 

REX  is  the  random  error  (meters)  associated  with 

fixing  the  East  coordinate  of  the  ascending 
radiosonde . 

REY  is  t>he  random  error  (meters)  associated  with 

fixing  the  North  coordinate  of  the  ascending 
radiosonde . 

RLE  is  the  random  error  (meters)  associated  with  the 
direct  measurement  of  the  launch  position. 

Further  problems  are  initiated  in  the  same  manner  as  described  for  RAWIN. 

Input  of  the  card  image  (?E0F  is  required  to  terminate  execution. 


3 . A  Example  Card  Input 


RAWIN: 


EXAMPLE  RAWIN  RUN  USING  15  ZONES  AND  3  ASCENT  RATES 

USE  INPUT  FILES  AR,  ZHHE ,  WF,  UV,  HRE,  AD712 

YES 

EXAMPLE  RAWIN  PROBLEM 

.03,  .05,  .03,  .05,  5. ,  0. ,  5. 

@E0F 
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RADAR: 


EXAMPLE  RADAR  RUN  USING  15  ZONES  AND  3  ASCENT  RATES 

USE  INPUT  FILES  AR,  ZHHE ,  WF,  UV,  HRE ,  AD712 

YES 

EXAMPLE  RADAR  PROBLEM 

.03,  .05,  .03,  .05,  0. ,  16. ,  5. ,  0. ,  5. 

@EOF 

NAVAID : 

EXAMPLE  NAVAID  RUN  FOR  15  ZONES  AND  3  ASCENT  RATES 

USE  INPUT  FILES  AR,  ZHHE,  UV,  WF 

YES 

EXAMPLE  NAVAID  PROBLEM 
100. ,  100.  ,  5. 

@E0F 
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4.0  FILE  INPUT 

Each  of  the  programs,  RAWIN,  RADAR,  and  NAVAID,  requires  data  to  be  inputted 
from  similar  sets  of  disk  files.  For  purposes  of  identification  each  required 
input  file  is  given  a  name.  Since  the  names  are  transparent  to  the  programs, 
the  user  may  rename  the  files  to  suit  his  own  needs. 

The  input  files  are  named,  respectively,  AR,  ZHHE,  UV,  WF,  KRE,  and  AD.  RAWIN 
and  RADAR  each  requires  all  of  these  files.  NAVAID  requires  only  AR,  ZHHE, 

UV,  and  WF.  In  this  section  the  files  are  described  in  the  same  order  that 
they  are  read  by  the  executing  program. 

The  logical  unit  which  each  program  uses  to  read  a  given  file  is  designated  by 
ION,  where  N  is  an  appropriate  integer.  For  example,  file  AR  is  read  from 
unit  102,  where  102  is  currently  assigned  the  value  2  If  the  user  desires, 
any  or  all  of  the  logical  units  may  be  reassigned  to  different  values  by 
changing  the  appropriate  data  statement  in  the  program. 

The  following  mnemonic  labels  are  used  here  to  describe  the  records  in  each 
file. 

LABEL  EXPLANATION 

IA  Balloon  ascent  rate  index.  IA  =  1,  NA 

NA  Total  number  of  balloon  ascent  rates.  NA  =  1,  2,  3,  or  4 

IZ  Ballistic  zone  index.  IZ  =  1,  NZ 

NZ  Total  number  of  ballistic  zones  in  the  complete  zone  structure.  For 

the  NATO  zone  structure  NZ  =  15.  However,  the  input  files  may  be 
constructed  for  any  value  of  NZ  between  1  and  30. 

IL  Ballistic  line  index.  IL  =  1,  NL 

NL  Total  number  of  ballistic  lines.  NL  =  NZ 

In  each  of  the  programs  various  arrays  involving  NZ  and/or  NL  are  currently 
dimensioned  to  allow  for  the  NATO  zone  structure  of  15  zones.  These  arrays 
must  be  redimens ioned  for  larger  zone  structures. 
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A.  1  File  AR 


File  AR  contains  the  balloon  ascent  rates. 

Required  by  RAWIN,  RADAR,  and  NAVAID 
Logical  unit  102  (currently  102  =  2) 

Total  number  of  records:  NA 
Number  of  data  items  per  record:  1 

RECORD  DATA  FORMAT 

IA  AR(IA)  Free  Field 

AR(IA)  is  the  ascent  rate  (meters/minute)  for  balloon 
IA. 

A. 2  File  ZHHE 

File  ZHHE  contains  zone  top  altitudes  and  the  bias  and  random  errors  asso¬ 
ciated  with  the  determination  of  these  altitudes. 

Required  by  RAWIN,  RADAR, and  NAVAID 
Logical  unit  103  (currently  103  =  3) 

Total  number  of  records:  NZ 

Required  number  of  data  items  per  record:  3 

RECORD  DATA  FORMAT 

IZ  Z(IZ),  BZ(IZ),  RZ(IZ)  Free  Field 

Z(IZ)  is  the  altitude  (meters)  at  the  top  of  zone 
IZ. 

BZ(IZ)  is  the  bias  error  (meters)  in  the  deter¬ 
mination  of  Z(IZ) . 

RZ(IZ)  is  the  random  error  (meters)  in  the  deter¬ 
mination  of  Z(IZ). 


4.3  File  UV 


File  U\  contains  the  zone  wind  protile 
Required  by  RAWIN,  RADAR,  and  NAVA  ID 
Logical  unit  1010  (currently  1010  =  10) 
Total  number  of  records:  NZ 
Required  number  of  data  items  per  record 


RECORD  DATA 

1Z  U(IZ),  V(IZ) 

U(IZ)  is  the  East  component  (meters/second)  of 

the  average  wind  velocity  in  zone  IZ. 

V(IZ)  is  the  North  component  (meters/second)  of 

the  average  wind  velocity  in  zone  IZ. 

4.4  File  WF 


File  WF  contains  the  zone  wind  weighting  factors. 
Required  by  RAWIN,  RADAR,  and  NAVAID 
Logical  unit  101  (currently  101  =  1) 

Total  number  of  records:  NL 

Required  number  of  data  items  per  record:  NZ 


RECORD  DATA 


IL  WF( 1 ) ,  WF(2) , . . . ,  WF(IZ) , . . . ,  WF(NZ) 


WF(IZ)  is  the  wind  weighting  factor  for  zone  IZ 
appropriate  to  ballistic  line  IL. 

Each  record  must  contain  NZ  values.  This  means  that 
for  IL  less  than  NL,  each  record  must  be  padded  at 
the  end  with  sufficient  zero  values  to  make  the  total 
number  of  data  items  in  the  record  equal  to  NZ. 


FORMAT 

Free  Field 


FORMAT 

Free  Field 
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4.5  File  HRE 


For  each  ascent  rate,  file  HRE  contains  values  of  positional  variables  relating 
to  the  launch  site  and  to  the  radiosonde  as  it  passes  each  zone  top. 

Required  by  RAWIN  and  RADAR 
Logical  unit  108  (currently  108  =  8) 

Total  number  of  records:  NA  *  (  NZ  +  1  ) 

Required  number  of  data  items  per  record:  2 

For  descriptive  purposes  we  divide  file  HRE  into  NA  successive  groups,  each 
group  corresponding  to  a  different  value  of  ascent  rate  index  IA  and  containing 
NZ  +  1  records.  The  following  is  a  description  of  records  within  any  one 
group. 

RECORD  DATA  FORMAT 

1  DISPL(IA),  AZL(IA)  (F10.0,  20X,  F10.2) 

DISPL(IA)  is  the  distance  (meters)  from  the  RAWIN  or 
radar  receiving  set  to  the  launch  site  for 
balloon  ascent  rate  index  IA. 

AZL(IA)  is  the  launch  azimuth  (degrees)  for  ascent 
rate  index  IA. 

1  +  IZ  E(IZ, IA) ,  A ( I Z , IA)  (20X ,  2F10.2) 

E(IZ,IA)  is  the  elevation  (degrees)  associated  with 
radiosonde  IA  at  the  top  of  zone  IZ. 

A  ( I Z , I A )  is  the  azimuth  (degrees)  associated  with 
radiosonde  IA  at  the  top  of  zone  IZ. 

For  each  ascent  rate  there  is  a  group  of  records  identical  to  the  group  listed 
above.  The  groups  must  be  appended  one  after  the  other  in  file  HRE  in  the 
same  order  that  the  ascent  rates  are  entered  in  file  AR.  There  are  no  blank 
records  or  other  markers  of  any  kind  between  the  groups. 
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The  X  fields  in  this  file  can  be  used  to  store  further  data  for  informational 
purposes.  For  example,  each  record  of  the  example  HRE  file  in  the  accom¬ 
panying  listing  contains  information  in  the  following  order:  horizontal 
distance,  slant  range,  elevation,  azimuth.  This  example  file  is  constructed 
for  3  ascent  rates  and  15  zones;  hence,  each  ascent  rate  group  contains  16 
records  to  yield  a  total  of  48  records  for  the  complete  file. 

4.6  File  AD 

File  AD  contains  bias  and  random  errors  in  elevation  associated  with  ground 
reflection.  These  errors  depend  on  a  number  of  factors,  including  the  antenna 
voltage  pattern  of  the  particular  receiving  set  under  consideration  and  the 
dielectric  constant  of  the  reflecting  surface.  The  example  AD  file  in  the 
accompanying  listing  is  mnemonically  called  AD712,  corresponding  to  a  seven 
foot  diameter  antenna  and  a  surface  dielectric  constant  of  12.  Program  LRDC 
may  be  used  to  generate  file  AD  in  the  appropriate  format. 

Required  by  RAWIN  and  RADAR 
Logical  unit  109  (currently  109  =  9) 

Total  number  of  records:  279 

Required  number  of  data  items  per  record  for  records  9  through  279:  2 

In  the  following  table  the  index  J  takes  on  the  values  J  =  1,  271. 

RECORD  DATA  FORMAT 

1--8  Miscellaneous  (///////) 

Records  1-8  contain  miscellaneous  information  des¬ 
cribing  the  parameters  which  were  used  to  create  the 
AD  file.  (See  example  listing.)  Program  LRDC  writes 
these  records  when  it  creates  this  file.  Although 
the  RAWIN  and  RADAR  programs  skip  these  records,  they 
must  be  present.  They  may  be  blank  if  the  user 
desires . 

8+J  BB(J) ,  RR( J)  (8X,  2F15.8) 
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BB(J)  is  the  bias  error  (degrees)  in  elevation 

due  to  ground  reflection  for  a  possible 

elevation  angle  of  FL0AT(J-l)/3  degrees. 

RR(J)  is  the  random  error  (degree)  in  elevation 

due  to  ground  reflection  for  a  possible 

elevation  angle  of  FL0AT(J-l)/3  degrees. 

In  other  words,  record  9  (J  -  1)  contains  error  data 
appropriate  to  0.0  degrees  elevation;  record  10 
(J  =  2)  contains  error  data  for  0.33  degrees  ele¬ 
vation;  record  11  (J  =  3)  contains  error  data  for 
0.67  degrees  elevation,  and  so  on  at  0.33  degree 
intervals.  Finally  record  279  (J  =  271)  contains 
error  data  for  90.0  degrees  elevation. 

If  desired  by  the  user,  the  X  field  in  each  record  may  contain  a  further  datum 

for  informational  purposes.  In  the  example  file  AD712  in  the  accompanying 

listing,  each  record  9  through  279  contains  information  in  the  following 
order:  possible  elevation  angle,  bias  error  due  to  reflection,  random  error 
due  to  reflection,  the  index  J.  (All  this  is  generated  by  program  LRDC.) 
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5 . 0  OUTPUT 


A  complete  run  of  RAWIN,  RADAR,  or  NAVAID  produces  two  output  print  files. 

The  program  writes  brief  output  to  logical  unit  106  (currently  106  =  6)  and 
more  extensive  output  to  logical  unit  1020  (currently  1020  =  20).  If  the  user 
desires,  the  logical  units  may  be  redesignated  to  different  values  by  changing 
the  appropriate  data  statement  in  each  program. 

5 . 1  Output  to  Unit  106 

In  demand  usage  this  output  is  directed  to  the  demand  terminal.  The  output  is 
formatted  for  printing  on  a  CRT  screen  which  can  display  80  characters  per 
line.  In  batch  mode  the  output  to  unit  106  is  redundant;  it  will  nevertheless 
be  directed  to  a  line  printer  unless  the  user  takes  steps  to  prevent  this. 

The  first  part  of  the  output  consists  of  solicitations  for  the  card  images 
described  in  Section  3  of  this  Users'  Manual. 

After  the  computations  for  any  given  problem  are  complete,  the  program  writes 
brief  results  to  unit  106.  These  results,  which  are  entirely  numerical,  may 
be  interpreted  from  the  following  table.  (In  this  table  the  term  'record'  is 
used  rather  loosely,  since  records  '2'  and  '3'  may  each  consist  of  several 
lines  of  output . ) 


RECORD 

DATA 

FORMAT 

1 

AR(  1 ) 

Free  Field 

AR(1)  is 
balloon  1. 

the 

ascent  rate  (meters/minute)  for 

'2' 

(CW(IL,  1) , 

IL=1, 

NL) 

(IX,  5F13.6) 

CW(IL,  1) 

is 

the  component  velocity  variance 

(knots**2)  in  the  ballistic  wind 
computed  for  ballistic  line  IL  and 
ascent  rate  index  1 . 
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'3' 


(SIGMA(IL) , 1 ) ,  IL=1 ,  ML) 


(IX,  5F13.6) 


SIGMA(IL,1)  is  the  standard  deviation  (knots) 
obtained  from  CW(IL,1J. 

For  a  given  problem  the  output  described  above  is 
repeated  for  each  balloon  ascent  rate  specified  in 
the  run. 

As  an  example,  consider  the  case  of  IS  ballistic  lines  and  one  ascent  rate. 

For  this  case  record  1  consists  of  one  printed  line  of  data  displaying  the 
balloon  ascent  rate.  Record  *2'  consists  of  three  printed  lines,  five  data 
items  per  line,  containing  the  values  of  the  CW  the  values  are  printed  from 
left  to  right  in  order  of  increasing  ballistic  line  number  IL.  Record  '3' 
also  consists  of  three  printed  lines,  five  data  items  per  line,  and  contains 
the  values  of  SIGMA;  the  values  are  printed  from  left  to  right  in  order  of 
increasing  IL.  Thus,  for  this  example,  the  total  number  of  printed  lines  of 
data  is  seven.  For  three  ascent  rates  21  lines  are  printed,  and  so  on. 

5 . 2  Output  to  Unit  1020 

In  either  demand  or  batch  mode,  this  output  is  meant  to  be  directed  to  a  line 
printer.  In  order  for  paging  and  line  spacing  commands  to  be  honored,  logical 
unit  1020  should  represent  an  alternate  print  file. 

If  the  user  entered  the  value  NO  on  input  card  3,  only  a  short,  self- 
explanatory  output  listing  is  produced.  It  consists  of  the  following: 

Comments  and  user  designated  file  names  appropriate  to  the  complete  run; 

Comments  and  input  data  for  each  individual  problem; 

2 

For  each  problem,  values  of  component  velocity  variance  (knots)  and 
standard  deviation  (knots)  computed  for  each  ballistic  line  for  each 
ascent  rate. 

If  the  user  entered  YES  on  card  3,  numerous  preliminary  and  intermediate 
results  are  also  outputted.  The  reader  is  referred  the  accompanying  listings 
for  examples  of  this  output  for  each  program. 


6.0  JOB  CONTROL 


The  job  control  run  streams  required  by  RAWIN,  RADAR,  and  NAVAID  are  similar 
for  all  three  programs.  The  main  difference  among  them  is  due  to  the  fact 
that  NAVAID  requires  fewer  input  files  than  the  other  two. 

Two  example  run  streams  for  RAWIN  are  shown  below  for  use  on  the  UN1VAC  1108. 
Many  variations  on  these  are  possible;  see  Reference  1.  Individual  variations 
for  RADAR  and  NAVAID  are  noted  at  the  end  of  each  run  stream. 


6.1  Initial  Batch  Run 


This  run  stream  achieves  the  following: 


Creates  and  lists  individual  temporary  input  data  files  from  cards,  where 
each  file  is  designated  by  an  appropriate  logical  unit  number; 

Creates  a  permanent  program  file  called  1NPUTFILES  and  copies  the  indi¬ 
vidual  temporary  data  files  into  separate  elements  in  INPUTFILES; 

Compiles  program  source  language  from  cards  and  creates  an  absolute 
executable  element; 

Places  symbolic,  relocatable,  and  absolute  elements  into  a  newly  created 
permanent  program  file  called  PROGRAM; 

Executes  the  program  for  a  complete  run  consisting  of  one  problem; 

Directs  all  output  to  a  specific  printed  labelled  PR3 . 


STEP  STATEMENT 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
I  1 
12 

13 

14 


@RUN  with  user  supplied 
@PASSWD  with  valid  user 
@SYM  PRINTS , 1 , PR3 
@ASG,CP  INP'JTFILES . 

@USE  I. .INPUTFILES. 
@ASG,CP  PROGRAM. 

@USE  P., PROGRAM. 

@ASG,CP  20. 

@ASG,T  2.  , / / / 32 
@ASG,T  3. , / / / 32 
@ASG,T  10.  ,///32 
@ASG ,T  1 .  ,  / / / 32 
@ASG , T  8. ,///32 
@ASG,T  9. ,///32 


options  and  information 
password 
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15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


@DATA,IL  2. 

--  AR  deck  goes  here  -- 
@END 

@DATA, IL  3. 

-*•  ZHHE  deck  goes  here  — 

@END 

@DATA,IL  10. 

--  UV  deck  goes  here  -- 
@END 

@DATA, IL  1. 

--  WF  deck  goes  here  -- 
@END 

@DATA,IL  8. 

--  HRE  deck  goes  here  -- 
@END 

@DATA, IL  9. 

--  AD  deck  goes  here  — 

@END 

@COPY , I  2. ,I.AR 
@C0PY,I  3. , I. ZHHE 
@COPY , I  10..I.UV 
@C0PY , I  1..I.WF 
@COPY,I  8. , I .HRE 
@COPY , I  9., I. AD 
@FREE  I. 

@FTN , IS  P.RAWIN 

--  RAWIN  source  deck  goes  here  -- 
@MAP , I  .P.RAWIN 
IN  P.RAWIN 
END 

@XQT  P.RAWIN 

Card  stating  user  comment  for  complete  RAWIN  run 
Card  designating  files  used  and/or  further  comment 
Card  requesting  or  rejecting  detailed  output 
Card  stating  user  comment  for  first  RAWIN  problem 
Card  containing  data  for  first  RAWIN  problem 
@E0F 

@FREE  20. 

@SYM  20. ,1,PR3 
@FIN 


Miscellaneous  Variations 


Steps  49  and  50  can  be  repeated  for  any  number  of  problems.  Insert 
additional  pairs  of  cards  after  step  50. 


The  L  option  of  each  @DATA  card  causes  a  listing  of  the  file  data  to  be 
produced.  If  this  listing  is  not  desired,  omit  the  L  option. 
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If  a  permanent  file  containing  data  elements  is  not  desired,  omit  Steps 
4,  5,  33--39 . 


If  a  permanent  file  PROGRAM  is  not  desired,  omit  Steps  6--7  and  replace 
Steps  40--45  with  the  following  sequence: 


@FTN,IS 

—  RAWIN  source  deck  goes  here  — 

@XQT 

If  only  file  creation  without  execution  is  desired,  omit  Steps  45--S3. 
Variations  for  RADAR 

The  job  control  sequence  for  RADAR  is  the  same  as  for  RAWIN.  In  the 
above  run  stream,  simply  replace  the  designation  RAWIN,  wherever  it 
occurs,  with  RADAR. 

Variations  for  NAVAID 

The  job  control  sequence  for  NAVAID  is  very  similar  to  that  for  RAWIN. 

In  the  above  run  stream,  simply  replace  the  designation  RAWIN,  wherever 
it  occurs,  wih  NAVAID.  Since  NAVAID  does  not  require  files  HRE  and  AD, 
the  following  steps  may  be  omitted:  13,  14,  27--32,  37,  38. 

b . 2  Typical  Demand  Run 

It  is  assumed  that  all  required  input  data  exist  in  appropriately  named  ele¬ 
ments  in  the  permanent  program  file  INPUTFILES.  Also,  it  is  assumed  that  an 
absolute  element  called  RAWIN  exists  in  file  PROGRAM.  Typically,  the  file 
PROGRAM  will  also  contain  symbolic  and  relocatable  elements. 

The  creation  of  appropriately  numbered  temporary  data  files  from  elements  in 
INPUTFILES  can  be  a  tedious  task  in  demand  mode.  To  facilitate  the  demand 
run,  it  is  assumed  that  the  user  has  previously  created  an  element  called,  for 
example,  RAWIN  in  a  program  file  called  ADDFILE.  The  contents  of  this  element 
are  shown  below. 


Contents  of  ADDFILE . RAWIN 


RECORD 

STATEMENT 

1 

@ASG,T  2.  ,///32 

2 

@ASG,T  3. ,///32 

3 

@ASG,T  10. , // / 32 

4 

@ASG,T  1.  ,///32 

5 

@ASG,T  8. , / / / 32 

6 

@ASG,T  9.  ,///32 

7 

@DATA, I  2. 

8 

@ADD , DP  I.AR 

.  UNIT  2 

9 

@END 

10 

@DATA , I  3. 

1 1 

@ ADD, DP  I.ZHHE 

.  UNIT  3 

12 

@END 

13 

@DATA , I  10. 

14 

@ ADD, DP  I.UV 

.  UNIT  10 

15 

@END 

16 

@DATA , I  1. 

17 

@ADD , DP  I.WF 

.  UNIT  1 

18 

@END 

19 

@DATA , I  8. 

20 

@ADD ,DP  I.HRE 

.  UNIT  8 

21 

@END 

22 

@DATA,I  9. 

23 

@ADD , DP  I. AD 

.  UNIT  9 

24 

@END 

An  element  ADDFILE . RADAR  would  be  identical  to  the  one  listed  above.  In  an 
element  ADDFILE .NAVAID  the  following  records  could  be  omitted:  5,  6,  19--24. 


Once  the  appropriate  ADDFILE  element  has  been  created,  it  can  be  used  with  any 
number  of  future  runs  in  demand  mode.  The  actual  demand  run  consists  of  the 
following  steps  (where  it  is  assumed  that  the  user  is  already  properly  signed 
on  to  the  terminal). 


STEP  STATEMENT 

1  @ASG,AZ  PROGRAM. 

2  @ASG, AZ  ADDFILE. 

3  @ASG , AZ  INPUTFILES. 

4  @USE  I. .INPUTFILES. 

5  @ASG ,CP  20. 

6  @ADD  ADDFILE. RAWIN 

7  @XQT  PROGRAM. RAWIN 

8  Entry  stating  user  comment  for  complete  RAWIN  run 

9  Entry  designating  input  files  and/or  further  comment 
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10  Entry  requesting  or  rejecting  detailed  output 

11  Entry  stating  user  comment  for  first  RAWIN  problem 

12  Entry  containing  data  for  first  RAWIN  problem 

1 3  @E0F 

14  (pFREE  20. 

15  @SYM  20. , 1 ,PR3 

16  @FIN  IF  DESIRED 

Miscellaneous  Variations 

Steps  11  and  12  may  be  repeated  in  sequence  for  any  number  of  problems. 

If  the  user  wants  to  test  changes  he  may  have  made  in  the  symbolic  element 
PROGRAM . RAWIN ,  he  can  easily  compile  this  element  and  execute  it  without 
retaining  new  relocatable  and  absolute  elements.  This  is  accomplished  by 
replacing  Step  7  with  the  following  sequence: 

@FTN , N  PROGRAM. RAWIN, TPF$. RAWIN 

@EGF 

@XQT 


The  complete  run  stream  listed  above  can  also  be  used  in  batch  mode  if 
the  following  three  cards  are  prefixed  to  it: 

@RUN  with  user  supplied  options  and  information 
(3PASSWD  with  valid  user  password 
@SYM  PRINTS, 1,PR3 

The  @SYM  card  listed  here  is  necessary  only  if  the  user  desires  to  direct 
the  PRINTS  output  (logical  unit  6)  to  the  specific  printer  PR3.  The  @FIN 
card  of  Step  16  is,  of  course,  required  in  batch  mode.  In  batch  usage, 
the  user  may  want  to  omit  Step  2  and  replace  Step  6  with  the  24  actual 
cards  listed  for  the  element  ADDFILE . RAWIN;  however,  this  is  not  necessary. 

Variations  for  RADAR 


The  job  control  language  for  RADAR  is  the  same  as  for  RAWIN.  In  the 
above  run  stream,  simply  replace  the  designation  RAWIN,  wherever  it 
exists,  with  RADAR.  In  Step  6,  the  element  ADDFILE . RAWIN  will  also  work 
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for  RADAR,  or  the  user  may  add  the  element  ADDFILE . RADAR  if  it  has  been 
created . 


Variations  for  NAVAID 

In  the  above  run  stream,  replace  the  designation  RAWIN,  wherever  it 
exists,  with  NAVAID.  It  is  assumed  that  an  element  ADDFILE. NAVAID  exists 
for  Step  6. 
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7.0  MODIFICATIONS 


Certain  easily  implemented  modifications  to  the  error  analysis  programs  are 
described  below: 

7 . 1  Logical  Units 

The  logical  units  numbers  used  in  input/output  operations  in  RAWIN,  RADAR,  and 
NAVAID  are  assigned  to  integer  variables  (beginning  with  characters  10)  in  a 
DATA  statement  in  each  program.  For  example,  in  RAWIN  we  have 

DATA  101,102,103,105,106,108,109,1010,1020/1,2,3,5,6,8,9,10,20/ 

The  user  may  change  the  assigned  values  in  the  DATA  statement  to  suit  his  con¬ 
venience  and/or  the  requirements  of  the  computer  system.  However,  it  is 
recommended  that  105  always  correspond  to  a  card  reader  or  terminal  keyin,  106 
to  a  line  printer  or  terminal  write,  and  1020  to  an  alternate  print  file. 

7 . 2  Redimensioning  of  Arrays 

Pertinent  arrays  in  RAWIN,  RADAR,  and  NAVAID  are  currently  dimensioned  to 
accomodate  15  zones ,  15  ballistic  lines,  and  four  balloon  ascent  rates.  The 
DIMENSION  statement  in  RAWIN,  for  example,  is  of  the  form: 

DIMENSION  AR(4) ,BB(271),RR(271) ,U(15) ,DISPL(4) ,AZL(4) ,V( 15) , W0( 15 ) 

1  ,WF(15)  ,7,(15)  ,A(  15,4)  ,CW(  15,4), DXDA(  15,4)  ,DXDE(15,4)  , 

2  DXDZ(15 ,4) ,DYDA( 15,4) ,DYDE(15,4) ,DYDZ( 15 ,4) , E( 15 , 4) , 

3  WV (15, 15), SI GMA (15,4), WX (15,4) ,VVY(15,4) ,W(15,15), 

4  WU(15, 15)  ,BZ(15),RZ(15)  ,'C0MM(20) 

Structures  of  15  zones  or  less  do  not  require  the  redimensioning  of  arrays. 

Of  course,  care  should  be  observed  in  preparing  the  various  input  data  files, 
as  described  in  Section  4.  (For  example,  each  of  the  programs  determines  the 
number  of  zones  and  ballistic  lines  from  the  number  of  records  in  file  WF.) 


Structures  of  up  to  30  zones  can  be  accommodated  by  the  programs.  For  struc¬ 
tures  containing  16  to  30  zones,  all  of  the  array  dimensions  currently  set  at 
15  must  be  changed  to  at  least  the  maximum  number  of  zones  in  the  structure. 
For  a  structure  of  26  zones,  for  example,  U(15)  becomes  U ( 26 1  ,  A ( 1 S , 4 )  becomes 
A(26,4),  W(15,15)  becomes  W(26,26),  and  so  on.  Array  dimensions  which  are  not 
currently  set  at  15  should  not  be  changed. 

7 . 3  Special  Changes  in  NAVAID 

If  the  user  desires,  any  of  the  values  assigned  in  the  following  DATA  state¬ 
ment  in  NAVAID  may  be  changed: 

DATA  TFIX,  TLM,  HINTL,  HINTH,  IZLO  / 1 . ,5 . ,200 . ,400 .  ,5/ 

where  the  variables  are  defined  below.  Note  the  following  restrictions: 

TFIX,  TLM,  HINTL,  and  HINTH  must  have  real  values  greater  than  zero;  IZLO  may 
be  zero  or  any  positive  integer. 


VARIABLE 

IZLO 

HINTL 

HINTH 


TFIX 

TLM 


EXPLANATION 

Highest  zone  for  which  height  interval  HINTL  is  used. 
Currently  IZLO  =  5 

Height  interval  (meters)  centered  at  zone  tops  and  uti¬ 
lized  in  fixing  the  East  and  North  coordinates  of  the 
ascending  balloon.  HINTL  is  used  for  IZ  =  1,  IZLO,  and 
HINTH  is  used  for  IZ  greater  than  IZLO.  Currently  HINTL 
=  200.  meters  and  HINTH  =  400  meters. 

Time  interval  (seconds)  between  successive  hyperbolic 
fixes  of  balloon  position.  Currently  TFIX  =  1.  second. 

Amount  of  time  (minutes)  used  in  fixing  the  initial  launch 
position.  Currently  TLM  =  5.  minutes. 
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8.0  PROGRAM  SEQUENCE  OF  OPERATIONS 


Each  program  performs  computations  and  input/output  operations  in  essentially 
the  same  sequence.  The  general  program  flow  is  given  below: 


STEP 

1 


2 

3 

4 

5 

6 


7 


8 


EXPLANATION 


Solicit  and  read  on  unit  105  three  card  images  containing  general 
documentation  and  instructions  for  the  complete  run.  Output  the 
general  documentation  to  unit  1020. 

Read  data  from  files  AR,  AHHE ,  and  UV  on  appropriate  logical  units. 

Read  zone  wind  weighting  factors,  line  by  line,  from  file  WF  on  unit 
101,  and  compute  the  weighting  arrays  W,  WU,  and  WV. 

Read  files  HRE  and  AD  on  the  appropriate  logical  units.  (Omit  this 
step  for  NAVAID. ) 

Compute  all  required  partial  derivatives  for  each  zone  for  each 
ascent  rate. 


Optionally  output  results  of  preliminary  computations  to  unit  1020. 
These  results  include  arrays  W,  WU,  and  WV,  as  well  as  arrays  con¬ 
taining  partial  derivatives. 


Solicit  and  read  on  unit  105  the  documentation  card  image  and  the 
data  card  image  for  specific  problem,  and  output  this  information  to 
unit  1020.  If  @E0F  is  read,  skip  to  Step  18. 


Do  through  Step  14  for  ascent  rate  index  IA  =  1 ,  NA. 


9  Compute  launch  component  errors  appropriate  to  IA. 

10  Do  through  Step  14  for  ballistic  line  IL  =  1 ,  NL. 

11  Compute  all  required  individual  error  sums  for  IL,  IA. 


12  Optionally  output  individual  error  sums  for  IL,  IA  to  unit  1020. 

IS  Compute  East  and  North  component  variances  in  ballistic  wind,  WX(IL, 

IA)  and  WY(IL,  IA),  respectively. 

14  Compute  the  component  velocity  variance  CW(IL,  IA)  in  ballistic 
wind  and  its  square  root  SIGMA(IL,  IA) . 

15  Optionally  output  to  unit  1020  the  arrays  WX  and  WY. 

Ic  Output  to  units  1020  and  106  the  arrays  CW  and  SIGMA. 
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17  Return  to  Step  17  for  next  problem. 

18  Terminate  execution. 


9.0  MNEMONICS 


In  all  three  error  analysis  programs  the  type  of  each  variable  is  in  general 
in  accordance  with  the  ASCII  FORTRAN  default  rule.  The  only  exceptions  to  the 
default  rule  involve  certain  explicitly  defined  character  type  variables 
containing  page  or  column  headings  which  are  written  to  output.  See  program 
listings . 

Numerous  variables  are  defined  elsewhere  in  this  Users'  Manual.  These  defini¬ 
tions  are  not  repeated  here.  For  convenience,  however,  reference  is  made 
below  to  the  sections  in  which  the  definitions  may  be  found.  Also,  the  various 
error  sums  and  required  partial  derivatives  are  discussed  under  separate 
headings.  Finally,  additional  miscellaneous  key  variables  are  defined  indi¬ 
vidual  ly . 

9 . 1  Variables  Defined  Elsewhere  In  This  Users'  Manual 

It  should  be  noted  that  the  physical  units  of  a  variable  may  change  during 
computation.  For  example,  the  program  reads  BA  in  degrees  and  later  converts 
the  units  to  radians. 

Variables  defined  in  Section  3: 

BA,  BEL,  BS,  COMM,  FE,  INTR ,  105,  RA,  RE,  REX,  REY,  RLE,  RLEA,  RLED,  RS. 
Variables  defined  in  Section  4: 

A(IZ,  IA),  AR(IA),  AZL(IA),  BB(J),  BZ(IZ),  DISPL(IA),  E(IZ,IA),  IA,  IL,  101, 
102,  103,  108,  109,  1010,  IZ,  NA,  NL,  NZ,  RR(J),  RZ(IZ),  U(IZ),  V(IZ),  WF(IZ), 
Z(IZ). 

Variables  defined  in  Section  5. 

ARU),  CW(IL,  1),  106,  1020,  SIGMA(IL,1). 
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Variables  defined  in  Section  7: 


HINTH,  HINTL,  IZLO,  TFIX,  TLM. 

9 . 2  Bias  and  Random  Error  Sums 

Current  values  of  the  individual  bias  and  random  error  sums  are  represented  by 
non  subscripted  variable  names  which  are  comprised  of  four  or  five  characters 
beginning  with  either  BE  or  RE.  See  Volume  I  of  this  report  of  the  definition 
of  each  of  these  sums. 

9 . 3  Partial  Derivatives  (RAWIN  and  RADAR) 

The  mnemonic  DpDq  represents  the  partial  derivative  of  p  with  respect  to  q. 
(However,  see  last  paragraph  under  the  current  subheading.)  In  the  Fortran 
code  p  and  q  are  written  as  characters  selected  from  the  following  list,  where 
all  distances  are  in  meters  and  all  angles  are  in  radians. 

p  OR  q  CHARACTER  EXPLANATION 

A  Azimuth  of  the  balloon-borne  radiosonde 

D  Distance  along  the  surface  of  the  earth  from  the  RAWIN  or  RADAR 

set  to  a  point  directly  below  the  ascending  radiosonde 

E  Elevation  of  the  radiosonde 

S  Slant  range  (RADAR  only) 

X  East  coordinate  of  radiosonde  position 

Y  North  coordinate  of  radiosonde  position 

Z  Altitude  of  the  radiosonde 

For  the  most  part,  values  of  the  partial  derivative  are  contained  in  arrays 
indexed  by  IZ  and  IA.  For  example,  DXDA(IZ,IA)  represents  the  partial  deri¬ 
vative  (meters/radian)  of  the  East  coordinate  of  radiosonde  position  with 
respect  to  azimuth,  appropriate  to  zone  IZ  and  ascent  rate  index  IA.  Non- 
subscripted  partials  represent  current  computational  values  only.  A  complete 
listing  of  partial  derivatives  utilized  in  RAWIN  and/or  RADAR  follows,  where 
IA  =  1,  NA,  and  IZ  =  1,  NZ;  partials  used  only  in  RADAR  are  marked  with  an 
a  s  te  r  i  s  k* : 
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DDDE ,  DDDZ ,  DXDA(IZ,IA),  DXDE(IZ,IA),  DXDZ(IZ,IA),  DYDA(IZ,IA),  DYDE(IZ.IA), 
DYDZ( IZ , IA) ,  DZDE(IZ.IA)*,  DXDS(IZ, IA)*. 


In  DO  loops  indexed  by  IA  and  IZ,  a  nonsubsrr ipteti  variable  of  the  form  DpDql 
is  used  as  a  convenience  to  hold  temporarily  the  array  value  DpDq(IZ.IA).  For 
example,  DXDAI  =  DXDA(IZ.IA). 

The  variables  DXDL,  DXDLA,  DYDL,  and  DYDLA  are  not  strictly  partial  derivatives. 
They  are  defined  below. 

9 . 4  Miscellaneous  Key  Variables 

Definitions  of  many  of  the  variables  used  in  RAWIN,  RADAR,  and  NAVAID  are 
given  elsewhere  in  this  Users’  Manual,  and  the  meanings  of  many  others  can  be 
readily  determined  from  an  examination  of  the  FORTRAN  code  itself.  Only 
certain  key  remaining  variables  are  defined  in  the  following  list.  Indication 
is  made  of  the  program(s)  in  which  each  variable  is  used. 


VARIABLE 


CVV ( I L , I A ) 


DXDLA 


DYDLA 


EXPLANATION 

Current  value  of  the  bias  error  (degrees,  radians)  in 
elevation  due  to  ground  reflection.  (RAWIN,  RADAR) 

2 

Component  velocity  variance  (knots)  in  ballistic  wind 
associtated  with  ballistic  line  IL  and  ascent  rate  IA. 
(RAWIN,  RADAR,  NAVAID) 

Current  value  of  the  product  (meters)  of  the  random  error 
in  launch  displacement  times  the  partial  derivative  of 
the  East  launch  coordinate  with  respect  to  launch  dis¬ 
placement.  (RAWIN,  RADAR) 

Current  value  of  the  product  (meters)  of  the  random  error 
in  launch  azimuth  times  the  partial  derivative  of  the  East 
launch  coordinate  with  respect  to  launch  azimuth.  (RAWIN, 
RADAR) 

Similar  to  DXDL  except  that  the  partial  derivative  of  the 
North  launch  coordinate  is  used.  (RAWIN,  RADAR) 

Similar  to  DXDLA  except  that  the  partial  derivative  of  the 
North  launch  coordinate  is  used.  (RAWIN,  RADAR) 


FNFL 


Conversion  factor,  degrees  to  radians.  (RAWIN,  RADAR) 

Real  variable  representing  the  number  of  hyperbolic  fixes 
of  balloon  launch  position.  (NAVAID) 


FNFZ 

G 

NFIXL 

NFIXZH 

NFIXZL 

Q 

R 

RER 

SIGMA(IL, IA) 

TOKNOT 

TORN  2 

W 

WX(IL,IA) 

WY ( IL ,  IA) 

W1 

W2 


Real  variable  representing  the  number  of  hyperbolic  fixes 
of  balloon  position  in  the  neighborhood  of  any  zone  top. 
(NAVAID) 

Current  value  of  the  angle  (radians)  subtended  at  the 
center  of  the  earth  by  the  measuring  set  and  the  ascending 
balloon.  (RAWIN,  RADAR) 


Number  of  hyperbolic  fixes  of  balloon  launch  position. 
(NAVAID) 


Number  of  hyperbolic  fixes  of  balloon  position  in  the 
neighborhood  of  zone  tops  greater  than  IZLO.  (NAVAID) 


Number  of  hyperbolic  fixes  of  balloon  position  in  the 
neighborhood  of  zone  tops  for  zones  IZ  =  1,  IZLO.  (NAVAID) 

Current  value  of  the  ratio  of  the  radius  of  the  earth  to 
the  sum  of  the  radius  of  the  earth  and  the  altitude  of  the 
ascending  balloon.  (RAWIN,  RADAR) 

Mean  radius  (meters)  of  the  earth.  (RAWIN,  RADAR) 


Current  value  of  the  random  error  (degrees,  radians)  in 
elevation  due  to  ground  reflection.  (RAWIN,  RADAR) 

Standard  deviation  (knots)  obtained  from  the  component 
velocity  variance  in  ballistic  wind  appropriate  ballistic 
line  IL  and  ascent  rate  IA.  (RAWIN,  RADAR,  NAVAID) 


Conversion  factor  by  division,  meters/second  to  knots. 
(RAWIN,  RADAR,  NAVAID) 

2 

Conversion  factor  by  division,  (meters/second)  to 
(knots).  (RAWIN,  RADAR,  NAVAID) 

Current-value  of  the  component  velocity  variance  (meters/ 
second;  in  ballistic  wind.  (RAWIN,  RADAR,  NAVAID) 

2 

Variance  (meters/second)  in  the  East  component  of  bal¬ 
listic  wind  for  ballistic  line  IL,  ascent  rate  IA. 

(RAWIN,  RADAR,  NAVAID) 

2 

Variance  (meters/second)  in  the  North  component  of  bal¬ 
listic  wind  for  ballistic  line  IL,  ascent  rate  IA.  (RAWIN, 
RADAR,  NAVAID) 

2 

Current  value  of  the  variance  (meters/second)  in  the  East 
component  of  ballistic  wind.  (RAWIN,  RADAR). 

2 

Current  value  of  the  variance  (meter/second)  in  the  North 
component  of  ballistic  wind. (RAWIN,  RADAR) 
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W(IZ,IL) 


WU( IZ , IL) 


WV(IZ,IL) 

WO(IL) 


A  defined  weighting  factor  (1/meter)  appropriate  to  bal¬ 
listic  line  IL  and  equal  to  the  wind  weighting  factor  per 
unit  zone  width  for  zone  IZ  minus  the  wind  weighting 
factor  per  unit  zone  width  for  zone  IZ  +  1  (RAWIN, 

RADAR,  NAVAID) 

A  defined  weighting  factor  (1/second)  appropriate  to  bal¬ 
listic  line  IL  and  equal  to  the  product  of  the  East  com¬ 
ponent  of  zone  wind  times  the  wind  weighting  factor  per 
unit  zone  width  (all  for  zone  IZ  +  1)  minus  the  product  of 
the  East  component  of  zone  wind  times  the  wind  weighting 
factor  per  unit  zone  width  (all  for  zone  IZ).  (RAWIN, 
RADAR,  NAVAID) 

A  defined  weighting  factor  (1/second)  similar  to  WU(IZ, 

IL)  except  that  the  North  component  of  zone  wind  is  used. 
(RAWIN,  RADAR,  NAVAID) 

A  defined  weighting  factor  (1/meter)  appropriate  to  bal¬ 
listic  line  IL  and  equal  to  the  negative  of  the  wind 
weighting  factor  per  unit  zone  width  for  zone  1.  (RAWIN, 
RADAR,  NAVAID) 
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10.0  PROGRAM  LRDC 


10.1  Introduction 

The  RAWIN  and  RADAR  models  require  as  input  an  estimate  of  the  error  in  ele¬ 
vation  angle  of  a  radiosonde's  position  due  to  the  fact  that  a  portion  of  the 
signal  reaching  the  RAWIN  or  RADAR  set  is  reflected  form  the  surface  of  the 
earth.  The  LRDC  model  provides  such  an  estimate. 

The  reader  is  referred  to  Reference  2,  for  example,  for  a  general  discussion 
of  tracking  by  the  sequential  lobing  (or  lobe  switching)  and  conical  scan 
techniques.  In  its  current  form  LRDC  is  a  simplified  modelling  of  the  sequen¬ 
tial  lobing  technique.  Reference  3  contains  a  useful  discussion  on  tracking 
errors  due  to  ground  reflections. 

The  intent  here  is  not  to  give  a  complete  description  of  LRDC,  which  is  still 
under  development,  but  rather  to  point  out  some  of  the  restrictions  associated 
with  its  use : 


The  antenna  patterns  for  the  upper  beam  and  lower  beam  switched  positions 
must  have  the  same  shape.  Further  development  will  be  required  to  remove 
this  restriction. 

In  order  to  find  the  attenuation  in  amplitude  of  a  reflected  (vertically 
polarized)  beam,  the  model  currently  assumes  that  the  reflectivity  is 
characterized  completely  by  the  angle  of  incidence  of  the  beam  at  the 
surface  of  the  reflecting  medium  and  be  the  real  dielectric  constant  of 
the  medium.  This  restriction  may  be  removed  in  further  development  to 
include  the  effect  of  signal  frequency  and  the  conductivity  of  the 
reflecting  medium. 


For  each  of  the  possible  elevation  angles  0.,  0.33,  0.67,  1.00,  1.33,  ..., 

90.0  degrees,  LRDC  computes  a  bias  error  in  elevation  and  a  random  error  in 
elevation  due  to  reflection.  The  computed  bias  error  associated  with  a  given 
elevation  is  typically  much  smaller  than  the  random  error  and  may  be  of  either 
sign. 


Program  LRDC  is  written  in  the  ASCII  FORTRAN  language.  It  may  be  executed  in 
either  demand  or  batch  mode. 
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10.2  Card  Input 


In  either  demand  or  batch  usage,  each  execution  of  program  LRDC  requires  one 
input  card  image  from  logical  unit  105  (currently  105  =  5). 


CARO 

DATA 

FORMAT 

1 

AS, 

DC,  SQA,  B 

Free  Field 

AS 

is  the  antenna  size  (dish  diameter)  in  any 
convenient  units.  This  parameter  is  used  solely 
for  documentation. 

DC 

is  the  dielectric  constant 
medium . 

of  the  reflecting 

SQA 

is  the  squint  angle  (degrees) 

appropriate  to  the 

antenna  pattern  of  the  tracking  device.  Typi¬ 
cally  the  upper  beam  antenna  voltage  pattern 
entered  from  file  AV  on  logical  unit  102  will 
already  have  the  squint  angle  built  into  it;  if 
this  is  the  case,  enter  0.  for  SQA. 

B  is  the  lower  beam  reduction  factor.  Enter  1.  if 
the  lower  beam  voltage  pattern  is  not  reduced  in 
magnitude  from  the  upper  beam  pattern. 


10.3  File  Input  ‘ 

Data  describing  the  upper  beam  antenna  voltage  pattern  are  inputted  from  file 

AV. 

Logical  unit  102  (currently  102  =  2). 

Total  number  of  records:  1080. 

Required  number  of  data  items  per  record:  4. 

In  the  following  table  I  takes  on  the  values  1=1,  1080. 

RECORD  DATA  FORMAT 


I  T 1 ( I ) ,  AV ( I ) ,  DAV(I) ,  DAV0(I)  Free  Field 

T 1  ( I )  is  the  angle  (degrees)  of  the  antenna 

beam  pattern  measured  from  the  forward 
antenna  axis.  It  must  have  the  value 
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given  by  T1(I)  =  FLOAT(  I  - 1  )/3 .  -180. 

In  other  words,  T 1 ( 1 )  =  -180.,  Tl(2) 
=  -179.67,  T1 (3)  =  -179.33,  .... 

T1 (541 )  =  0.00,  ....  11(1080)  = 

179.67,  all  in  degrees. 

AV(I)  is  the  absolute  value  of  the  upper 

beam  voltage  corresponding  to  T 1 ( I )  . 
Note  that  the  values  of  the  AV(I)  are 
normalized  such  that  the  maximum  value 
of  the  set  of  all  the  A V  is  equal  to 
1.  (This  maximum  value  will  typically 
be  located  near,  but  not  necessarily 
precisely  at,  index  I  =  541.) 

DAV(I)  is  an  antenna  voltage  difference  given 
by  DAV(I)  =  AV(I+1)  -AV(I). 

DAV0(I)  is  an  antenna  voltage  difference  given 
by  DAV0(I)  =  (AV(I+1)  -AV(I-l) )/2 . 


The  variables  AV(I),  DAV(I),  and  DAVO(I)  should  be  entered  to  as  many  signi¬ 
ficant  figures  as  is  feasible. 


10.4  Output 

Program  LRDC  writes  output  to  logical  units  103  and  106.  Currently  103  =  3, 
and  106  =  6.  These  units  can  be  redesignated  to  different  values  by  changing 
the  appropriate  data  statement  in  LRDC. 

10.4.1  Output  to  Unit  103 

The  output  to  unit  103  can  be  used  to  create  an  AD  file  suitable  for  use  by 
programs  RAW1N  and  RADAR. 

Total  number  of  records:  279. 

Number  of  data  items  per  record  for  records  9  through  279:  5. 

In  the  following  table  J  assumes  the  values  J  =  1,  271. 

RECORD  DATA  FORMAT 

1--8  Miscellaneous  Various 

These  records  repeat  the  card  input  data,  dis¬ 
play  the  computed  elevation  offset  (degrees), 
and  write  column  headings.  (See  AD712  listing.) 
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8+J 


T1(J),  ELAV(J) ,  SIGMA(J) ,  RMS(J),  J 


(1X,F7.2,3F15.8,I5) 


T1(J)  is  the  possible  target  elevation  angle 

(degrees)  given  by  T1(J)  =  FLOAT(J- 

l)/3.  (The  program  redefines  the  array 
T1  from  its  input  values  given  in 
section  10.3  to  the  values  noted  here 
for  output . ) 

ELAV(J)  Is  the  computed  mean  elevation  error 
(degrees)  due  to  reflection  for  angle 
T1(J). 

SIGMA(J)  is  the  computed  standard  deviation 

(degrees)  of  elevation  errors  due  to 
reflection  for  angle  T 1 ( J ) . 

RMS(J)  Is  the  computed  root  mean  square  error 
(degrees)  in  elevation  due  to  reflec¬ 
tion  for  angle  Tl(J). 

J  Is  the  counting  index  J  =  1,271. 

The  values  ELAV(J)  and  SIGMA(J)  are  taken,  respectively,  to  be  the  bias  and 
random  error  due  to  reflection  for  elevation  angle  T1(J).  The  quantities 
RMS(J)  and  J  are  outputted  for  informational  purposes  only. 


The  output  to  unit  103  is  in  the  correct  from  of  an  AD  file  for  input  to 
programs  RAWIN  and  RADAR.  See  the  example  file  AD712  in  the  accompanying 
1 isting. 


10.4.2  Output  to  Unit  106 

In  batch  usage  the  output  to  unit  106  will  result  in  a  printed  listing.  In 
demand  mode  this  output  is  directed  to  the  demand  terminal. 


Output  to  unit  106  is  similar  to  the  ''uput  to  unit  103.  However,  the  following 
differences  should  be  noted:  The  first  eight  records  in  the  above  table  are 

omitted  ;  the  counting  index  J  is  replaced  by  an  iteration  index  appropriate 
to  each  elevation  angle. 
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10.5  Job  Control 


Shown  below  are  two  example  LRDC  run  streams  for  use  on  the  UN  I  VAC  1108. 
Numerous  variations  on  these  are  possible.  See  Reference  1. 

10.5.1  Initial  Batch  Run 


This  run  stream  achieves  the  following: 


Creates  from  cards  a  data  element  AV7  in  a  previously  catalogued  file  called 
INPUTFILES; 

Compiles  LRDC  source  language  form  cards  and  creates  an  absolute  executable 
element ; 

Places  symbolic,  relocatable,  and  absolute  elements  into  a  previously  cata¬ 
logued  file  called  PROGRAM; 

Executes  LRDC  for  the  case  of  dielectric  constant  of  10.; 

Places  the  resulting  output  to  unit  3  into  the  element  INPUTFI LES . AD7 10 . 

STEP  STATEMENT 

1  @RUN  with  user  supplied  options  and  information 

2  (9PASSWD  with  valid  user  password 

3  @SYM  PRINTS , 1 , PR3 

4  @ASG . A  INPUTFILES. 

5  @USE  I.,  INPUTFILES. 

b  (pASG ,  A  PROGRAM. 

7  <8USE  P .  ,  PROGRAM . 

8  @ASG,  T  2. 

0  @ASG,  T  3. 

1 0  @DATA ,12. 

11  --  AV  7  deck  goes  here  -- 

12  @ END 

13  @C0PY,  I  2. ,  I . AV7 

14  @FTN ,  I S  P .  LRDC 

15  --  LRDC  deck  goes  here  -- 

lb  @MAP , I  , P . LRDC 

17  IN  P. LRDC 

18  END 

19  @XQT  P . LRDC 

20  7 . ,  10. ,  0.  ,  1 . 

21  @ELT , ID  I.AD710 

22  @ADD,D  3. 

23  @END 

24  @FIN 
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It  is  assumed  that  the  input  data  from  file  AV7  exists  in  an  element  called 
INPUTFILES . AV7 .  It  is  also  assumed  that  an  absolute  executable  element  called 
PROGRAM. LRDC  exists. 


The  program  is  executed  for  an  input  dielectric  constant  of  12.  the  resulting 
ouput  to  logical  unit  3  is  catalogued  in  a  data  file  called  AD712  and  is  also 
placed  in  an  element  INPUTFILES . AD7 12 . 

After  the  user  is  properly  signed  on  to  the  terminal,  one  possible  sequence  of 
commands  is  as  follows: 


STEP  STATEMENT 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 
1G 


@ASG , AZ  INPUTFILES. 
@USE  I. .INPUTFILES. 
@ASG,AZ  PROGRAM. 
@USE  P.,  PROGRAM. 
@ASG , CP  AD712. 

(RISE  3. ,AD7 12 . 
@ASG,T  2. 

(RATA,  I  2. 

@ADD,D  I.AV7 
(SEND 

@XQT  P . LRDC 
7.  ,  12.  ,  0.  ,  1. 

@ELT , ID  I.AD712 
0ADD.D  3. 

@END 

@FIN.  IF  DESIRED 


REFERENCES 


Computer  Branch,  NROD,  WSMR  UNIVAC  1108  L'ser  Guide,  Technical  Report  No. 
69. 


Skolnik,  Merrill  I.,  1962,  Introduction  To  Radar  Systems,  McGraw-Hill, 
NY.  . ~ 


Barr,  William  C. ,  and  Peterson,  Arnold  C.,  1977,  "Wind  Measuring  Accuracy 
Test  of  Meteorological  Systems",  ECOM-5831,  Atmospheric  Sciences  Labora¬ 
tory,  US  Army  Electronics  Command,  White  Sands  Missile  Range,  NM. 


COMPUTER  LISTINGS 


The  following  pages  contain  computer  listings  of  programs  RAWIN',  RADAR,  NAVAID, 
and  LRDC .  Example  input  tiles  are  also  listed  as  well  as  output  from  sample 
runs.  (Due  to  its  length  an  example  File  AV  is  not  listed.) 
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